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We elaborate on the consequences of the factorization of the transfer matrix of any lossless mul- 
tilayer in terms of three basic matrices of simple interpretation. By considering the bilinear trans- 
formation that this transfer matrix induces in the complex plane, we introduce the concept of 
multilayer transfer function and study its properties in the unit disk. In this geometrical setting, 
our factorization translates into three actions that can be viewed as the basic pieces for understand- 
ing the multilayer behavior. Additionally, we introduce a simple trace criterion that allows us to 
classify multilayers in three types with properties closely related to one (and only one) of these 
three basic matrices. We apply this approach to analyze some practical examples that are typical 
representatives of these types of matrices. 

PACS numbers: 



I. INTRODUCTION 

Layered media play an important role in many appli- 
cations in modern optics, especially in relation to optical 
filters and the like. Therefore, it is not surprising that 
the topics covered in most of the textbooks on the subject 
use a mixture of design, manufacture, and applications; 
dealing only with the basic physics needed to carry out 
practical computations [jjj . 

However, for a variety of reasons, layered media have 
physical relevance on their own ^ [| . As any linear sys- 
tem with two input and two output channels, any multi- 
layer can be described in terms of a 2 x 2 transfer matrix. 
In fact, it has been recently established that for a lossless 
multilayer this transfer matrix is an element of the group 
SU(1,1) j|, H). From this perspective, it is precisely the 
abstract composition law of SU(1,1) the ultimate respon- 
sible for the curious composition law of the reflection and 
transmission coefficients || [?J . 

This purely algebraic result is certainly remarkable. 
But, as soon as one realizes that SU(1T) is also the ba- 
sic group of the hyperbolic geometry it is tempting 
to look for an enriching geometrical interpretation of the 
multilayer action. Moreover, given the role played by ge- 
ometrical ideas in all branches of physics, particularly in 
special relativity, it is easy to convince oneself that this 
approach might provide deeper insights into the behavior 
of a multilayer in a wider unifying framework that can 
put forward fruitful analogies with other physical phe- 
nomena. 

Accordingly, we have proposed || to view the action 
of any lossless multilayer as a bilinear transformation 
on the unit disk, obtained by stereographic projection 
of the unit hyperboloid of SU(1,1). This kind of bilin- 
ear representations have been discussed in detail for the 
Poincare sphere in polarization optics jl0[ [ll] , for Gaus- 
sian beam propagation [ |l2j , and are also useful in laser 
mode- locking and optical pulse transmission fll3[. 



In spite of these achievements, the action of an ar- 
bitrary lossless stack could still become cumbersome to 
interpret in physical terms. In fact, in practice it is usual 
to work directly with the numerical values of a matrix 
obtained from the experiment, which cannot be directly 
related to the inner multilayer structure. To remedy this 
situation, we have resorted recently juj to the Iwasawa 
decomposition, which provides a remarkable factoriza- 
tion of the matrix representing any multilayer (no mat- 
ter how complicated it could be) as the product of three 
matrices of simple interpretation. 



At the geometrical level, such a decomposition trans- 
lates directly into the classification of three basic actions 
in the unit disk, which are studied in this paper, that 
can be considered as the basic bricks for understanding 
multilayers. Moreover, we have also shown that the 
trace of the transfer matrix allows for a classification of 
multilayers in three different types with properties very 
close to those appearing in the Iwasawa decomposition. 



In this paper we go one step further and exploit this 
new classification to study several practical examples 
that are representatives of each type. This shows the 
power of the method and, at the same time, allows for 
a deeper understanding of layered media. As a direct 
application, we treat the outstanding case of symmetric 
multilayers, finding a precise criterion for zero-reflectance 
conditions. Nevertheless, we stress that the benefit of 
this formulation lies not in any inherent advantage in 
terms of efficiency in solving problems in layered struc- 
tures. Rather, we expect that the formalism presented 
here could provide a general and unifying tool to an- 
alyze multilayer performance in an elegant and concise 
way that, additionally, is closely related to other fields of 
physics, which seems to be more than a curiosity. 
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II. TRANSFER MATRIX FOR A LOSSLESS 
MULTILAYER 



We first briefly summarize the essential ingredients of 
multilayer optics we shall need for our purposes |l0| . 
The configuration is a stratified structure, illustrated in 
Fig. 1, that consists of a stack of 1, . . . , j, . . . , m, plane- 
parallel lossless layers sandwiched between two semi- 
infinite ambient (a) and substrate (s) media, which we 
shall assume to be identical, since this is the common 
experimental case. Hereafter all the media are supposed 
to be lossless, linear, homogeneous, and isotropic. 

We consider an incident monochromatic linearly polar- 
ized plane wave from the ambient, which makes an angle 
9q with the normal to the first interface and has ampli- 
tude Ea . The electric field is either in the plane of 
incidence (p polarization) or perpendicular to the plane 
of incidence (s polarization). We consider as well another 
plane wave of the same frequency and polarization, and 
with amplitude E^ ' , incident from the substrate at the 
same angle 9q. @ 

As a result of multiple reflections in all the interfaces, 
we have a backward-traveling plane wave in the ambient, 
denoted Ei~\ and a forward-traveling plane wave in the 
substrate, denoted E^ + K If we take the field amplitudes 



as a vector of the form 



E 



E (-) 



(1) 



which applies to both ambient and substrate media, then 
the amplitudes at each side of the multilayer are related 
by the 2x2 complex matrix M as , we shall call the mul- 
tilayer transfer matrix, in the form 

E Q = M as E s . (2) 

The matrix M as can be shown to be of the form B 



Ras/T as 1/T* S 





a 


(3 ' 




(3* 


* 

a 



(3) 



where the complex numbers 

Ras = \Ras\ exp(zp), 



T as = \T as \exp(iT), (4) 



are, respectively, the overall reflection and transmission 
coefficients for a wave incident from the ambient. Note 
that we have det M as = +1, which is equivalent to 
|-Ras| 2 + \T a s\ 2 — 1, and then the set of lossless multilayer 
matrices reduces to the group SU(1,1), whose elements 
depend on three independent real parameters. 

The identity matrix corresponds to T as = 1 and R as — 
0, so it represents an antireflection system without trans- 
mission phase shift. The matrix that describe the over- 
all system obtained by putting two multilayers together 
is the product of the matrices representing each one of 
them, taken in the appropriate order. So, two multilay- 
ers, which are inverse, when composed give an antireflec- 
tion system. 



If we denote by R sa and T sa the overall reflection and 
transmission coefficients for a wave incident from the sub- 
strate (physically this corresponds to the same multilayer 
taken in the reverse order) one can check that 



T as Tsa ~ RasRsa = exp(z2-r), 



Rsa = -#Q S ex P(«2T), 



(5) 



which is a generalization of the well-known Stokes rela- 
tions H for the overall stack. 

On the other hand, it is worth noting that, while the 
picture of a multilayer taken in the reverse order is clear, 
at first sight it is not so easy to imagine the inverse of 
that multilayer. However, using Eqs. (|5|), one can obtain 
that 

M^ 1 = MJ a , (6) 
which remedies this drawback. 



III. A BASIC FACTORIZATION FOR 
MULTILAYERS: THE IWASAWA 
DECOMPOSITION 

Many matrix factorizations have been considered in 
the literature Jl8], ||(| , the goal of all of them being to 
decompose a matrix as a unique product of other matri- 
ces of simpler interpretation. Particularly, given the es- 
sential role played by the Iwasawa decomposition, both in 
fundamental studies and in applications to several fields 
(especially in optics), one is tempted to investigate also 
its role in multilayer optics. 

Without embarking us in mathematical subtleties, the 
Iwasawa decomposition is established as follows pl[ : any 
element of a (noncompact semi-simple) Lie group can be 
written as an ordered product of three elements, taken 
one each from a maximal compact subgroup K, a maxi- 
mal Abelian subgroup A, and a maximal nilpotent sub- 
group N. Furthermore, such a decomposition is global 
and unique. 

For a lossless multila yer matrix M QS 6 SU(1,1), the 
decomposition reads as [pL4j] 



where 



K(0) 
A(0 
N(iz) 



„ = K(0)A(0N(i/) , 

exp(#/2) 

exp(-#/2) 

cosh(£/2) i sinh(£/2) 
-i sinh(£/2) cosh(£/2) 

l-iv/2 v/2 
u/2 l + iu/2 



(7) 



(8) 



The parameters <f>, £, and v are given in terms of the 
elements of the multilayer matrix by 

(j)/2 = arg(a + i/3) , 
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HI = ln(l/|a + »/3|) , 
v/2 = Re{a[3*)/\a + if3\ 2 



(9) of zero R as by the simple condition 



and their ranges are and — 2tt < <fi < 2tt. There- 

fore, given a priori limits on a and fj (i.e., on T as and 
R as ), one could easily establish the corresponding limits 
on (f>, £ and v , and vice versa. 

Now, we can interpret the physical action of the matri- 
ces appearing in Eq. (Q). K(0) represents the free prop- 
agation of the fields E in the ambient medium through 
an optical phase thickness of (f>/2, which reduces to a 
mere shift of the origin of phases. Alternatively, this can 
be seen as an antireflection system. The second matrix 
A(£) represents a symmetric multilayer with real trans- 
mission coefficient TX = sech(£/2) and reflection phase 
shift pa = ±7r/2. There are many ways to get this perfor- 
mance, perhaps the simplest one is a Fabry-Perot system 
composed by two identical plates separated by a trans- 
parent spacer. By adjusting the refractive indices and the 
thicknesses of the media one can always get the desired 
values (see Section 7). Finally, the third matrix, N(i/), 
represents a system having Tn = cos(tn) exp(iT>j) and 
Rfi = sin(ri\i) exp(ir|\|), with tan(ri\i) = v/2. The sim- 
plest way to accomplish this task is by an asymmetrical 
two-layer system. 



IV. A REMARKABLE CASE: SYMMETRIC 
MULTILAYERS 

So far, we have discussed various properties of arbi- 
trary lossless multilayers. However, the particular case 
of symmetric structures has deserved a lot of attention. 
In this Section, we wish to show how the Iwasawa de- 
composition provides a nice tool to deal with them. 

We recall that for a symmetric stack the reflection and 
transmission coefficients are the same whether light is 
incident on one side or on the opposite side of the mul- 
tilayer. In consequence, one has R as = R sa and the 
generalized Stokes relations (^|) give then the well-known 



result 1 22, 



p-T = ±n/2. (10) 

This implies that the element (3 in the transfer matrix (|^) 
is a pure imaginary number. Therefore, the matrix M as 
depends only on two real parameters, which translates 
into the fact that 0, f, and v [see Eqs. @] are not in- 
dependent. In fact, a straightforward calculation shows 
that they must fulfill the constraint 



v = (e s - l)tan(0/2); 
and thus R as can be written as 

= ^ tanh(g/2)[tan(0/2)-i] 
1 -ztanh(£/2)tan(0/2)' 



Ras 



(11) 



(12) 



Particular care has been paid to the characterization 
of zero reflectance conditions for these symmetrical sys- 
tems p5|, p6l . Equation dT5) allows us to express the loeii 



C = 0, 



(13) 



which, by Eq. ([fl]) implies v = and imposes that in this 
case M as reduces trivially to a matrix K (</!>). 

The stability of these nonreflecting configurations has 
been studied by Lekner pr| . Indeed, by using a conti- 
nuity argument, he has shown that "almost all partial 
reflectors with symmetric profiles which are close in pa- 
rameter space to a profile which has reflectivity zeros, 
will also have reflectivity zeros" . We intend to show how 
our formalism allows for a more precise criterion. 

To this end, let us assume a symmetric multilayer sat- 
isfying initially the condition (|l^). Now, suppose that 
some parameter(s) (refractive index, thickness, angle of 
incidence, ...), we shall generically denote by £, is var- 
ied. Obviously, admissible variations must preserve the 
symmetry of the system. 

The variation of t induces changes in R as , and so in <p 
and £. The new multilayer will have also zero R as if the 
parameters satisfy dR as /d£ = 0; that is, 



dR as 
d£ 



dR n 



00 



dR n 



d( 



di 



d£ 



(14) 



Using Eq. (|12|), one gets that dR as /d(j)\^ = o is identically 
zero, while dR as /d^\^ = o never vanishes. We conclude 
then that the condition we are looking for is 



*=0. 

dl 



(15) 



This result fully characterizes the partial reflectors in- 
voked by Lekner, and can be of practical importance for 
the design of robust antireflection systems. 



V. MULTILAYER TRANSFER FUNCTION IN 
THE UNIT DISK 

In many instances (e.g., in polarization optics [To|), we 
are interested in the transformation properties of quo- 
tients of variables rather than on the variables them- 
selves. In consequence, it seems natural to consider the 
complex numbers 



£(+)' 



(16) 



for both ambient and substrate. From a geometrical 
viewpoint, Eq. (^|) defines a transformation of the com- 
plex plane C, mapping the point z s into the point z a , 
according to 



*[M QS ,z s 



(3* +a*z s 
a + (3z s 



(17) 



Thus, the action of the multilayer can be seen as a func- 
tion z a = f{z s ) that can be appropriately called the mul- 
tilayer transfer function The action of the inverse 



matrix 
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These bilinear transformations define an action of the 
group SU(1,1) on the complex plane C. The complex 
plane appears then decomposed in three regions that re- 
main invariant under the action of the group: the unit 
disk, its boundary and the external region. p7[ 

The Iwasawa decomposition has an immediate trans- 
lation in this geometrical framework, and one is led to 
treat separately the action of each one of the matrices 
appearing in this decomposition. To this end, it is worth 
noting that the group SU(1,1) that we are considering ap- 
pears always as a group of transformations of the complex 
plane. The concept of orbit is especially appropriate for 
obtaining an intuitive picture of the corresponding ac- 
tion. We recall that, given a point z, its orbit is the 
set of points z' obtained from z by the action of all the 
elements of the group. In Fig. 2 we have plotted some 
typical orbits for each one of the subgroups of matrices 
K (</>), A(£), and IM(V). For matrices K(0) the orbits are 
circumferences centered at the origin and passing by z. 
For A(£), they are arcs of circumference going from the 
point +i to the point — i through z. Finally, for the ma- 
trices N (v) the orbits are circumferences passing through 
the point +i and joining the points z and — z* . 



VI. TRACE CRITERION FOR 
CLASSIFICATION OF MULTILAYERS 

To go beyond this geometrical picture of multilay- 
ers, let us introduce the following classification: a ma- 
trix is of type K when [Tr(M as )] 2 < 4, is of type A 
when [Tr(M as )] 2 > 4, and finally is of type N when 
[Tr(M as )] 2 = 4. To gain insight into this classification, 
let us also introduce the fixed points |2f| of a transfer ma- 
trix as the points in the complex plane that are invariant 
under the action of M as ; i.e., 



$[M as ,z], 



whose solutions are 



-ilm(a) ± y/[Re(a)] 2 - 1 




Since we have 



Tr(M as ) = 2Re(a) 



2 COST 



(18) 



(19) 



(20) 



one can easily check that the matrices of type K have 
two fixed points, one inside and other outside the unit 
disk, both related by an inversion; the matrices of type 
A have two fixed points both on the boundary of the unit 
disk and, finally, the matrices of type N have only one 
(double) fixed point on the boundary of the unit disk. 

Now the origin of the notation for these types of ma- 
trices should be clear: if one consider the Iwasawa de- 
composition (^), one can see that the matrices K(0) are 
of type K with the origin as the fixed point in the unit 
disk, matrices A(£) are of type A with fixed points +i 



and —i and matrices N(^) are of type N with the double 
fixed point +i. Of course, this is in agreement with the 
orbits shown in Fig. 2. 

To proceed further let us note that by taking the con- 
jugate of M QS with any matrix C G SU(1,1) we obtain 
another multilayer matrix; i.e., 



M, 



C 



(21) 



such that Tr(M as ) = Tr(M as ). The fixed points of M as 
are then the image by C of the fixed points of M as . If we 
write the matrix C as 



C 



Cl 



C2 
C? 



(22) 



the matrix elements of M as (denoted by carets) and those 
of M as are related by 



a|ci| 2 -a*|c 2 | 2 -2Tm(/3cic;), 



(i = /3c 2 -/3*c^-2icic 2 lm(a). 



(23) 



For our classification viewpoint it is essential to remark 
that if a multilayer has a transfer matrix of type K, A, or 
N, one can always find a family of matrices C such that 
M as in Eq. (||) is just a matrix K((f>), A(£), or N(i/), re- 
spectively. The explicit construction of this family of ma- 
trices is easy: it suffices to impose that C transforms the 
fixed points of M os into the corresponding fixed points of 
K (</>), A(£), or N(^). By way of example, let us consider 
the case when M os is of type K and its fixed point inside 
the unit disk is Zf. Then, one should have 

$[CM as C-\ 0] = $[CM as , Zf ] = $[C, z f ] = 0. (24) 

Solving this equation one gets directly 



Cl 



1 



: exp(iS), 



C2 



(25) 



where 5 is a real free parameter. The same procedure 
applies to the other two cases. 

Since the matrix M as belongs to one of the subgroups 
K(</>), A(£), or N(is) of the Iwasawa decomposition, and 
all these subgroups are, in our special case, Abelian and 
uniparametric, we have that 



M/ii)M as (^2) = M as (/ii +/i 2 ), 



(26) 



where /i represent the adequate parameter 0, £, or v. 
Therefore, when dealing with a periodic layered system, 
whose matrix is obtained as the iVth power of the basic 
period, we have 



M 



-lMiV 



f s (/i)C = C- 1 M as (iV M )C, (27) 



where M as is now the matrix of the basic period. This 
is a remarkable result E9] , and our procedure highlights 
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that it does not depend on the explicit form of the basic 
period. 

We wish to point out that the trace criterion has been 
previously introduced |3(| Q to treat the light propa- 
gation in periodic structures. In fact, in that approach 
the values of the trace separate the band-stop from the 
band-pass regions of a period stratification. 

Finally, in order to broaden the physical picture of this 
classification, let us transform Eq. (|J) by the unitary ma- 
trix 



y/2 



1 i 
i 1 



Then, we can rewrite it alternatively as 



(28) 



(29) 



where the new field vectors £ and the new multilayer 
matrix M a s are obtained by conjugation by U. 

One can easily check that det.M as = +1 and all its 
elements are real numbers. Therefore, M a s belongs to 
the group SL(2,R) that underlies the structure of the 
celebrated ABCD law in first-order optics |3^, |33 , 34 35) . 

By transforming by U the Iwasawa decomposition (m) , 
we get the corresponding one for SL(2,R), which has been 
previously worked out |36fl : 



M as = lC(<l>)A{Z)Af{v), 



(30) 



where 



/C(0) = 
A(0 = 



cos(^/2) sin(0/2) 
- sin(0/2) cos(0/2) 

ex P (£/2) 
e X p(-£/2) 

1 0" 

v 1 



(31) 



The physical action of these matrices is clear. Let us 
consider all of them as ABCD matrices in geometrical 
optics that apply to position x and momentum p (direc- 
tion) coordinates of a ray in a transverse plane. These 
are the natural phase-space variables of ray optics. Then 
/C(0) would represent a rotation in these variables, A(£) 
a magnifier that scales x up to the factor m = exp(£/2) 
and p down by the same factor, and M{v) the action of 
a lens of power v (34j| . 

In the multilayer picture, £^ can be seen as the corre- 
sponding x, while can be seen as the corresponding 
p. Then, the key result of this discussion is that when 
the multilayer transfer matrix has [Tr(M as )] 2 lesser than, 
greater than or equal to 4 one can find in a direct way 
a family of matrices that gives a new vector basis such 
that the action of the multilayer, when viewed in such 
a basis, is exclusively rotationlike, or magnifierlike, or 
lenslike [fill. 



VII. SIMPLE EXAMPLES AND CONCLUDING 
REMARKS 

It seems pertinent to conclude by showing how our ap- 
proach works in some practical examples. Perhaps, the 
best way of starting is to consider the simplest layered 
structure one can imagine: a single film sandwiched be- 
tween the same ambient and substrate media. In spite 
of its simplicity, it contains the essential physical ingre- 
dients of multilayer optics. 

In consequence, we consider a single transparent film 
(medium 1) of refractive index m and thickness d\ em- 
bedded in air (medium 0). For this system we have g§ 



T as — Tow — 



(1 -rg 1 )cxp(-i/3i) 
1 - r§! exp(-z2/?i) : 



(32) 



where r i is the Fresnel reflection coefficient at the in- 
terface 01 and /?i = (27midi cos#i)/A is the plate phase 
thickness (here A is the wavelength in vacuo of the inci- 
dent light and 0\ is the refraction angle). Accordingly, 
we get 



[Tr(Moio)] 2 = 4cos 2 /3i < 4, 



(33) 



and the equality holds only in resonance conditions; i.e., 
when |Tbio| = | cos /3i | = 1. In consequence, the matrix 
of a single film is always of type K . 

Let us consider now two films (1 and 2) described by 
the matrices Moio and M020, respectively. The compound 
system obtained by putting them together is described 
by the product of these matrices M ioM 2o [@, H 0- In 
conclusion, since any layered stack can be viewed as the 
composition of single films, this shows that any multilayer 
matrix is generated by the product of matrices of type 
K. Take into account that the product of two matrices 
of type K (or A or N) can have trace lesser than, greater 
than or equal to 4. 

On the other hand, as we have stated at the end of 
Section 3, to get a pure matrix A(£) one should consider 
a Fabry-Perot-like system formed by two identical plates 
(each one of them with phase thickness (3\ ) separated by 
a spacer of air with phase thickness /?2 ■ If we take as ini- 
tial condition that in the substrate z s = 0.4exp(— in/ 3), 
then a standard calculation gives T as and R as , and from 
them we obtain the value z a — —0.44 + 0.49 i, with the 
parameters indicated in Fig. 3. Obviously, from these 
(experimental) data alone we cannot infer at all the pos- 
sible path for this discrete transformation. 

However, the Iwasawa decomposition remedies this se- 
rious drawback: from the geometrical meaning discussed 
before, and once we know the values of <fi, £, and v [that 
are easily computed from Eqs. (0)] we get, by the or- 
dered application of the matrices K(<f>), A(£), and N(j/), 
that the trajectory from z s to z a is well defined through 
the corresponding orbits, as shown in Fig. 3. 

Moreover, and this is the important moral we wish to 
extract from this simple example, if in some experiment 
the values of z s and z a are measured, one can find, no 
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matter how complicated the multilayer is, in a unique 
way, the three arcs of orbits that connect the initial and 
final points in the unit disk. 

In Fig. 4 we have plotted the values of the parameters 
cj>, £, and v for this system when /?2 is varied between 
and 7r. We have also plotted the values of [Tr(M as )] 2 . It 
is evident that the system can be of every type depending 
on the value of @2 ■ The marked points determine special 
behaviors in agreement with Eq. ([ll]) for symmetrical 
systems: for the left one, (j> = v = and the system is 
represented by a pure matrix A(£); in the right one, £ = 
v = and then it is represented by a matrix K(0); i.e., it 
is antireflection stack with T as = exp(— i<fi/2). Note that 
this system can never be represented by a matrix N(z^), 
because it is symmetric. 

To show the characteristic properties of an asymmetric 
system, we consider, as indicated in Section 3, the sim- 
plest one constituted by a two-layer stack made of a glass 



plate (with phase thickness j3i) coated with a film of zinc 
sulphide (with phase thickness /3 2 ). In Fig. 5 we have 
plotted the values of 4>, £, and v, as well as [Tr(M as )] 2 , 
when (3\ is fixed and fii is varied between and it. From 
our previous analysis it is clear that only in the marked 
point we have = ^ = and [Tr(M QS )] 2 = 4, so the 
system is represented by a pure matrix N(^). Contrary 
to the previous example, this system can never be repre- 
sented by a matrix A(£), because it is asymmetric. 

In summary, we expect that the geometrical scenario 
presented here could provide an appropriate tool for an- 
alyzing and classifying multilayer performance in an ele- 
gant and concise way that, additionally, could be closely 
related to other fields of physics. 

Luis L. Sanchez-Soto's e-mail address is 
lsanchez@eucmax.sim.ucm.es 
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FIG. 1: Wave vectors of the input \E^ and Ei '] and output 



(_) and E { a +) ] fields in a multilayer sandwiched between two 



FIG. 2: Plot of several orbits in the unit disk of the elements 
of the Iwasawa decomposition K(0), A(£), and N(^) for the 
group of multilayer transfer matrices. 



identical semi-infinite ambient and substrate media. 



FIG. 3: Geometrical representation in the unit disk of the 
action of a symmetric system made up of two identical plates 
(m = 1.7, di — 1 mm, 8o = 7r/4, A = 0.6888 fim and s- 
polarized light) separated by a spacer of phase thickness /?2 = 
3 rad. The point z s — 0.4 exp(— J7t/3) is transformed by the 
system into the point z a = —0.44 + 0.49i. We indicate the 
three orbits given by the Iwasawa decomposition and, as a 
thick line, the trajectory associated to the multilayer action. 
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FIG. 4: Plot of the values of [Tr(M as )] 2 and of the parameters 
<j>, £, and v in the Iwasawa decomposition for the same system 
as in Fig. 3, as a function of p2- 



FIG. 5: Same plot as in Fig. 4 but for an asymmetric system 
made up of a glass plate (m = 1.5 and /3i = 2.75 rad) coated 
with a zinc sulphide film (ri2 = 2.3). 
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